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The ALICE collaboration has made the first measurement at the LHC of J/ y/ photoproduction in 
ultra-peripheral Pb-Pb-collisions at v / 5nn= 2.76 TeV. The J/i// is identified via its dimuon decay in 
the forward rapidity region with the muon spectrometer for events where the hadronic activity is 
required to be minimal. The analysis is based on an event sample corresponding to an integrated 
luminosity of about 55 /ib~ 1 . The cross section for coherent J/y production in the rapidity interval - 
3.6 <y < -2.6 is measured to be dcr^/dy — 1.00±0.18(stat)+o2g(syst) mb. The result is compared 
to theoretical models for coherent production and found to be in good agreement with those 
models which include nuclear gluon shadowing. 
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Two-photon and photonuclear interactions at unprecedentedly high energies can be studied in ultra- 
peripheral heavy-ion collisions (UPC) at the LHC. In such collisions the nuclei are separated by impact 
parameters larger than the sum of their radii and therefore hadronic interactions are strongly suppressed. 
The cross sections for photon induced reactions remain large because the strong electromagnetic field of 
the nucleus enhances the intensity of the virtual photon flux, which grows as Z 2 , where Z is the charge 
of the nucleus. The virtuality of the photons is restricted by the nuclear form factor to be of the order 
l/Rm 30 MeV/c (R is the radius of the nucleus). The physics of ultra-peripheral collisions is reviewed 
in (DEI- 

Exclusive photoproduction of vector mesons, where a vector meson but no other particles are pro- 
duced in the event, is of particular interest. Exclusive production of J/y in photon-proton interactions, 
7 + p— >J/y + p, has been successfully modelled in perturbative QCD in terms of the exchange of two 
gluons with no net-colour transfer Q- Experimental data on this process from HERA have been used 
to constrain the proton gluon-distribution at low Bjorken-x (3). Exclusive vector meson production in 
heavy-ion interactions is expected to probe the nuclear gluon-distribution Q, for which there is con- 
siderable uncertainty in the low-* region [6). A J/y produced at rapidity y is sensitive to the gluon 
distribution at x = (M T/ y/ y^nn) exp(±;y) at hard scales Q 2 pa 0. The two-fold ambiguity in x is 

due to the fact that either nucleus can serve as photon emitter or photon target. At the forward rapidities 
studied here (-3.6 < y < -2.6), the relevant values of x are 10~ 2 and « 10~ 5 , respectively. 

Exclusive p° j8l and J/y (9j production have been studied in Au-Au collisions at RHIC. The p° is too 
light to provide a hard scale, and the J/y analysis suffered from very low statistics, so no conclusions 
concerning nuclear shadowing were made from these studies. Exclusive J/y production has also been 
studied by the CDF collaboration in proton-antiproton collisions at the Tevatron iflOl . The availability 
of such measurements has led to an increase in interest in ultra-peripheral collisions, stimulating several 
new model calculations. 

In this Letter, the first LHC results on exclusive photoproduction of J/y vector mesons are presented. 
J/y mesons produced in Pb-Pb collisions at ^Ann = 2.76 TeV have been measured at forward rapidities 
through their dimuon decay. Exclusive photoproduction can be either coherent, where the photon cou- 
ples coherently to all nucleons, or incoherent, where the photon couples to a single nucleon. Coherent 
production is characterized by low vector meson transverse momentum ((pr) — 60 MeV/c) and the tar- 
get nucleus normally does not break up. Incoherent production, corresponding to quasi-elastic scattering 
off a single nucleon, is characterized by a somewhat higher transverse momentum ((px) — 500 MeV/c) 
and the target nucleus normally breaks up, but except for single nucleons or nuclear fragments in the 
very forward region no other particles are produced. This analysis is focussed on coherently produced 
J/y mesons. The experimental definition of coherent production, which must take into consideration 
also the finite detector resolution, is here pj < 0.3 GeV/c. The measured cross section is compared to 
model predictions El H21 E3 El • 

The ALICE detector consists of a central barrel placed inside a large solenoid magnet (B = 0.5 T), 
covering the pseudorapidity region |tj| < 0.9 lfT31 . and a muon spectrometer covering the range -4.0< 
tj <-2.5. The spectrometer consists of a ten interaction length (A/) thick absorber filtering the muons, 
in front of five tracking stations containing two planes of cathode pad multi-wire proportional chambers 
(MWPC) each, with the third station placed inside a dipole magnet with a fBdl = 3 Tm integrated 
field. The forward muon spectrometer includes a triggering system, used to select muon candidates 
with a transverse momentum larger than a given programmable threshold. It has four planes of resistive 
plate chambers (RPC) downstream of a 1.2 m thick iron wall (7.2 A/), which absorbs secondary punch- 
through hadrons from the front absorber and low momentum muons from % and K weak decays. This 
analysis uses the VZERO counters for triggering and event selection. These consist of two arrays of 32 
scintillator tiles each, covering the range 2.8< r\ <5.1 (VZERO-A, on the opposite side of the muon 
arm) and -3.1 <rj <-1.7 (VZERO-C) and positioned at z = 329 cm and z = -87 cm from the interaction 
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point, respectively. Finally, two sets of hadronic Zero-Degree Calorimeters (ZDCs) are located at 116 
m on either side of the Interaction Point. These detect neutrons emitted in the very forward region, for 
example neutrons emitted following electromagnetic dissociation ifToll . 

The analysis presented in this publication is based on a sample of events collected during the 2011 
Pb-Pb run, selected with a special trigger (FUPC) set up to select UPC events in which a particle is pro- 
duced at forward rapidities. The integrated luminosity corresponds to about 55 [lb~ y ; the instantaneous 
luminosity during the Pb-Pb run varied from 10 24 cm _2 s _1 , with 1 interacting bunch per beam at the 
beginning of the data taking period to 2.5 • 10 26 cm~ 2 s _1 , with 336 interacting bunches per beam at the 
end of the data taking period. 

The purpose of the FUPC trigger is to select events containing two muons in an otherwise empty detec- 
tor, from two-photon production (77 — > /i + jit~) or from J/y decay, and it requires the following event 
characteristics: 

(i) a single muon trigger above a 1 GeV/c /^-threshold; 

(ii) at least one hit in the VZERO-C detector; 

(iii) no hits in the VZERO-A detector. 

A total of 3.16 x 10 6 events were selected by the FUPC trigger. 

The offline event selection used in a previous J/y analysis ifTTl was modified to account for the typical 
experimental signatures of ultra-peripheral processes, i.e. only two tracks in the spectrometer and very 
low J/ y transverse momentum. The following selection criteria were applied: 

(i) two reconstructed tracks in the muon arm; 

(ii) owing to the multiple scattering in the front absorber, the DCA (distance between the vertex and the 
track extrapolated to the vertex transverse plane) distribution of the tracks coming from the interaction 
vertex can be described by a Gaussian function, whose width depends on the absorber material and is 
proportional to 1 / p, where p is the muon momentum. The beam induced background does not follow 
this trend, and was rejected by applying a cut on the product p x DCA, at 6 times the standard deviation 
of the dispersion due to multiple scattering and detector resolution. The additional dispersion due to the 
uncertainty on the vertex position (not measurable in UPC events) is negligible in comparison and does 
not affect the value of the cut; 

(iii) at least one of the muon track candidates were required to match a trigger track above the 1 GeV/c 
/^-threshold in the spectrometer trigger chambers; 

(iv) both tracks pseudorapidities within the range -3.7< T]i 2 <-2.5, to match the VZERO-C acceptance; 

(v) the tracks exit from the absorber in the range 17.5 cm< /? a t, s <89.5 cm, delimiting the two homo- 
geneous parts of the absorber covering the angular acceptance of the spectrometer (/? a bs is the radial 
coordinate of the track at the end of the front absorber); 

(vi) dimuon rapidity to be in the range -3.6< y <-2.6, which ensured that the edges of the spectrometer 
acceptance were avoided; 

(vii) two tracks with opposite charges; 

(viii) only events with a neutron ZDC signal below 6 TeV on each side were kept. In the present data 
sample, this cut does not remove any events with a J/y produced with a transverse momentum below 
0.3 GeV/c, but reduces hadronic contamination at higher py, 

(ix) dimuons to have pj < 0.3 GeV/c and invariant mass 2.8 < M- mv < 3.4 GeV/c 2 ; 

(x) VZERO offline timing compatible with crossing beams. 

After applying these selections, 117 J/y candidates remained. The effect of the cuts on the statistics is 
listed in Table [T] 

The acceptance and efficiency of J/i/A-reconstruction were calculated using a large sample of coherent 
and incoherent J/y events generated by STARLIGHT lfT8l and folded with the detector Monte Carlo 
simulation. STARLIGHT simulates photonuclear and two-photon interactions at hadron colliders. The 
simulations for exclusive vector meson production and two-photon interactions are based on the models 
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Table 1: Summary of the applied data cuts (see text). 



Selection 



Number of remaining events 



Triggered events 
Two charged tracks 
p x DCA cut 

RPC muon chamber matching 

-3.7<T]i j2 < -2.5 

Absorber thickness (R a bs) 

-3.6 <y < -2.6 

1 OS dimuon 

Neutron ZDC < 6 TeV 

p T < 0.3 GeV/c and 2.8 < M inv < 3.4 GeV/c 2 

VZERO offline timing 



3,161,675 
432,422 
26,958 
10,712 



5,100 
5,095 
4,919 
3,209 
817 
122 
117 



in ifTTI and |[T9l , respectively. These models are discussed further below. 

The residual misalignment and the time-dependent conditions of the tracking and trigger chamber com- 
ponents were taken into account in these simulations. The trigger chamber efficiencies were computed 
from the data and used in the global efficiency calculation. A separate simulation was performed for 
each run, in order to take into account the slight variations in run conditions during the data taking. The 
product of the acceptance and efficiency corrections (Acc x £)t/^ was calculated as the ratio of the num- 
ber of the simulated events that satisfy the event selection in Table 1 to the number of generated events 
within -3.6 < y < -2.6. The final values for the combined acceptance and efficiency were found to be 
16.6% and 14.3% for coherent and incoherent J/y, respectively. The relative systematic error coming 
from the uncertainties on the trigger chamber efficiencies used in these simulations amounts to 4%. In 
addition, the muon reconstruction efficiency has been evaluated both in data and simulations, in a way 
similar to that described in IfTTI . and a 6% relative systematic uncertainty on the (Acc x £)j/y corrections 
was assigned to account for the observed differences. 

In order to evaluate the systematic error on the acceptance coming from the generator choice, the ac- 
ceptance was computed from a parameterization of the results on coherent J/y production in 0. It 
was also calculated by modifying the rapidity distribution in STARLIGHT and letting it vary between 
a flat distribution and a distribution consistent with the model with the steepest slope (AB-MSTW08, 
see below for definition) over the range -3.6 < y < -2.6. The differences in acceptance between the 
methods were below 3%, which was taken into account in the systematic error calculation. It is assumed 
in these calculations that the J/y is transversely polarized. Transverse polarization is expected for a 
quasi-real photon from s-channel helicity conservation. This has been confirmed experimentally for ex- 
clusive J/y production in y+p — > J/y + p interactions lf20ll2Tl and for exclusive p° photoproduction in 
heavy-ion collisions (H. Owing to the low pj of the J/y, the calculations are insensitive to the choice of 
reference frame (here the helicity frame was used), and the polarization axis effectively coincides with 
the beam axis. 

Activity in the central barrel was checked for events with invariant mass in the range 2.8 < M- mv < 
3.4 GeV/c 2 . No events with more than one tracklet in the Si-Pixel (SPD) detector were found. The 
events with one tracklet (6 out of 117) were not removed, as this level of activity is consistent with the 
expected background from random combinations of noise hits. 

The invariant mass distribution for opposite sign (OS) muon pairs with 2.2 < M- mv < 4.6 GeV/c 2 is 
shown in Fig. 1. A J/y peak is clearly visible in the spectrum, on top of a continuum coming from 
77 — > jU + jU~. Only two like-sign dimuon pairs are in the invariant mass range 2.2 < M mw < 4.6 GeV/c 2 , 
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at 2.3 GeV/c 2 and 2.8 GeV/c 2 . The combinatorial background is therefore estimated to be <2% at 90% 
confidence level in the invariant mass range 2.8 < M mv < 3.4 GeV/c 2 . 

The J/ y yield was obtained by fitting the dimuon invariant mass spectrum in the range 2.2 < M im < 4.6 
GeV/c 2 with an exponential function to describe the underlying continuum, and a Crystal Ball function 
E2l to extract the J/y signal. The Crystal Ball tail parameters ((Xcb and n) were fixed to values obtained 
from simulations. The y' was not considered in the fitting procedure as its contribution is negligible in 
the dimuon channel. The central mass value from the fit is 3.123 ±0.011 GeV/c 2 , which is within 2.4a 
(0.8%) of the known value of the J/y mass and compatible with the absolute calibration accuracy of 
the muon spectrometer. The width, 84±14 MeV/c 2 , is in agreement with the Monte Carlo simulations. 
The extracted number of J/ys is A^eid = 96 =L 12(stat) ± 6(syst). The systematic error on the yield 
(6%) was obtained by varying the Crystal Ball tail parameters. The exponential slope parameter of the 
continuum is -1.4±0.2 GeV _1 c 2 in good agreement with the corresponding Monte Carlo expectation 
(-1.39±0.01 GeV _1 c 2 ). This is an additional indication that there is no unexpected background in the 
invariant mass and pj range considered. 

The fraction of the J/y mesons coming from the decay of y — > J/y + anything was estimated by 
simulating a sample of coherently produced y s with STARLIGHT, using PYTHIA ll23l to simulate 
their decay into J/y. The detector response was simulated as described above. The contribution from 
incoherently produced y is expected to give a negligible contribution for pj < 0.3 GeV/c and was not 
considered. The y polarization and its relative transfer to the J/y affect the muon angular distribu- 
tion and thus their acceptance. Unlike the directly produced J/y discussed above, the polarization of 
J/ys coming from y decays cannot easily be predicted, since the polarization of the original y can be 
shared between the J/y and the other daughters in different ways. The y decay was therefore simu- 
lated by assuming the following J/ y polarizations: (i) no polarization (NP); (ii) full transverse (T), and 
(iii) full longitudinal (L). The J/y fraction coming from y decay for a given polarization P, f£ can be 
written as: 



P _ <V BR (V ^ W + anything) • (Acc x e)^ J/y 
D Oj /v , • (Acc x e) J/¥ 

where the (Acc x e)j/y, and (Acc x were computed for pj < 0.3 GeV/c. 

According to STARLIGHT, the ratio between the y and J/y coherent photoproduction cross sections 
is 0.19 giving f^ p = 11.9%, f£ = 9.3%, = 16.8%. The cross sections ratio is significantly lower in 
the pQCD inspired model Q, 0.087. This changes the above fraction, giving f$ p = 5.5%, f£ = 4.3%, 
= 7.9%. The estimates for fjj thus range from 4.3% to 16.8%. The best estimate was taken as the 
middle of this range with the extremes providing the lower and upper limits, giving /b=(l 1±6)%. 

The dimuon pj distribution integrated over 2.8 < M- mv < 3.4 GeV/c 2 is presented in Fig. 2. The clear 
peak at low pj is mainly due to coherent interactions, while the tail extending out to 0.8 GeV/c comes 
from incoherent production. In addition, the high-px region may still contain a few hadronic events, 
which makes it difficult to extract the incoherent photoproduction cross section from these data. To 
estimate the fraction (/}) of incoherent over coherent events in the region pj < 0.3 GeV/c, the ratio 
cw/Ocoh> weighted by the detector acceptance and efficiency for the two processes, was calculated, giv- 
ing fi = 0.12 when ai nc /a co h was taken from STARLIGHT, and fj = 0.08 when the model in Q was 
used. Four different functions were used to describe the pj spectrum: coherent and incoherent photo- 
production of J/y, J/y from y decay, and two-photon production of continuum pairs. The shapes for 
the fitting functions (Monte Carlo templates) were provided by STARLIGHT events folded with the de- 
tector simulation. The relative normalization was left free for coherent and incoherent photoproduction. 
The contribution from the y was constrained from estimate above (/b=( 1 1 ±6)%), and the two-photon 
contribution was determined from the fit of the continuum in Fig. 1 . In the fit, the incoherent process 
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is constrained mainly in the region 0.5 < pj < 0.8 GeV/c, where the other three processes are negli- 
gible. As this pt region (not used in the J/y signal extraction) is likely to suffer from some hadronic 
background, the fit can only provide an upper limit on //. The result is // = 0.26 ±0.05, about a factor 
2 larger than the estimate from the theoretical models quoted above. We conclude by taking the middle 
value of the three calculations as the best estimate of fj, and the other two results as lower and upper 
limits, respectively, giving fj = 0.121q^. 

Finally, the total number of coherent J/ys is calculated from the yield extracted from the fit to the 
invariant mass distribution by 

resulting in Nf?* = 78 ± 10(stat) % (syst). 

The coherent J / y differential cross section is given by: 



N coh 

-2U. (3 ) 

dy (Accxe) ]/v -e t n g -BR(J/Y — > n+n^-Jf^-Ay ' 

where NjJ^ is the number of J/y candidates from Eq. |2l (Acc x e)j/y, corresponds to the acceptance 
and efficiency of the muon spectrometer, as discussed above, and e t rig is the VZERO trigger efficiency. 
BR(J/y — ► /i + jU _ ) = 5.93% is the branching ratio for J/y decay into muons GH, Ay = 1 the rapidity 
interval bin size, and ££- mX the total integrated luminosity. During the 2011 Pb-Pb run the VZERO de- 
tector was optimised for the selection of hadronic Pb-Pb collisions, with a threshold corresponding to an 
energy deposit above that from a single minimum ionizing particle (MIP). The distribution of the signal 
produced by a MIP crossing the 2 cm thick VZERO scintillator has a Landau shape. To get an accurate 
simulation of the efficiency for low multiplicity events with this threshold setting, would require an al- 
most perfect reproduction of the Landau by the MC simulation. Therefore we used the QED continuum 
pair production for normalization. 

In addition to exclusive J/y, the FUPC trigger selected 77 — )■ A t+ A t events, which are very similar 
to coherent J/y decays in terms of kinematics and associated event characteristics. This reaction is a 
standard QED process, which in principle can be calculated with high accuracy. The fact that the photon 
coupling to the nuclei is Z^fa (with Z = 82 here) rather than just ^fa increases the uncertainty of the 
contribution from higher order terms. Predictions exist where this effect is negligible |[25l . However, 
other studies obtained a 16% reduction in the cross section from higher order terms in Pb-Pb collisions 
at the LHC l26l . There is also an uncertainty associated with the minimum momentum transfer and 
the nuclear form factor |[27l . Two-photon production of /x + /i _ -pairs from STARLIGHT was used to 
determine the trigger efficiency |fl9l . The cross sections from STARLIGHT for two-photon production 
of e + e~ and jU + /l~ pairs have previously been compared with results from STAR [28] and PHENIX 
respectively. The predictions from STARLIGHT have been found to be in good agreement with the 
experimental results. These results, however, have uncertainties of about 25 to 30%. In the absence of 
high precision measurements constraining the model, and taking into account the outstanding theoretical 
issues mentioned above, the uncertainty in the STARLIGHT two-photon cross section is estimated to be 
20%. 

The cross section for 77 — > pL + pL~ can be written in a similar way to Eq. 3 and the ratio of the two is 
independent of luminosity and of the trigger efficiency: 



J/y 



da 



coh 

j/y 



1 



coh 



J/y 



dy BR(J/y ^ jx+fx- 



rr 



(Acc x e 
(Acc x e 



rr 
Vy 



'rr 



Ay 



(4) 
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Table 2: Summary of the contributions to the systematic uncertainty for the integrated J/ y cross section mea- 
surement. The error for the coherent signal extraction includes the systematic error in the fit of the invariant mass 
spectrum and the systematic errors on fr> and //, as described in the text. 



Source 


Value 


Theoretical uncertainty in a 7 y 


20% 


Coherent signal extraction 


-14 /0 


Reconstruction efficiency 


6% 


RPC trigger efficiency 


5% 


J/y acceptance calculation 


3% 


Two-photon e + e~ background 


2% 


Branching ratio 


1% 


Total 


+24 % 

-26 /0 



where Nyy was obtained by counting the number of events in the invariant mass intervals 2.2 < Mj nv < 
2.6 GeV/c 2 (N rY = 43) and 3.5 < M inv < 6 GeV/c 2 (N 77 = 15), to avoid contamination from the J/y peak. 
To determine Oyy STARLIGHT pT9l was used. The cross section for dimuon invariant mass between 2.2 
< M inv < 2.6 GeV/c 2 or 3.5 < M im < 6 GeV/c 2 , dimuon rapidity in the interval -3.6 < y < -2.6, and 
each muon satisfying -3.7<t/i 2< -2.5 is a yr = 17.4 /lb (a yr = 13.7 jttb and Oy 7 = 3.7 /lb for the low and 
high invariant mass intervals, respectively). The (Acc x e)yy for events satisfying the same selection was 
calculated using events from STARLIGHT folded with the detector simulation as described above. The 
data cuts applied to the Monte Carlo sample were the same as those applied for the J/y data analysis, 
resulting in a (Acc x e) yy of 42.1% (37.9% for 2.2 < M im < 2.6 GeV/c 2 and 57.5% for 3.5 < M inv < 
6 GeV/c 2 ). 

A possible source of inefficiency comes from correlated QED pair production, i.e. interactions which 
produce both a J/i|/ and a low mass e + e~-pair (the latter has a very large cross section), with one of the 
electrons hitting the VZERO-A detector and thus vetoing the event. This effect was studied with data, in 
a sample collected with comparable luminosity by a control trigger, requiring a coincidence of at least 
two muons in the muon arm trigger with hits in both the VZERO-A and VZERO-C. Two J/y events 
were found in this sample, giving an upper limit on the inefficiency smaller than 2%. 

Since the kinematic distributions of the muons from J/y decays and 77 processes are different, the 
systematic uncertainties on the corresponding (Acc x e) corrections coming from the uncertainties on 
the muon trigger and reconstruction efficiencies do not exactly cancel out in equation [4] In order to 
account for this effect, a 50% correlation factor has been estimated, conservatively, when comput- 
ing the systematic uncertainty on the ratio. The sources of the systematic error are summarized in 
Table |2] The final result is a differential cross section for coherent J/y production of da^/dy = 

1.00±0.18(stat)l^(syst) mb. 

The cross section is compared with calculations from various models (51 [TT] [I2j [131 in Fig. 3. The 
five models calculate the photon spectrum in impact parameter space in order to exclude interactions 
where the nuclei interact hadronically. The differences between the models come mainly from the way 
the photonuclear interaction is treated. The predictions can be divided into three categories: 

i) those that include no nuclear effects (AB-MSTW08, see below for definition). In this approach, all 
nucleons contribute to the scattering, and the forward scattering differential cross section, da/dt at t = 
(t is the momentum transfer from the target nucleus squared), scales with the number of nucleons squared, 
A 2 ; 

ii) models that use a Glauber approach to calculate the number of nucleons contributing to the scattering 
(STARLIGHT, GM, and CSS). The reduction in the calculated cross section depends on the total J/y- 
nucleon cross section; 
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iii) partonic models, where the cross section is proportional to the nuclear gluon distribution squared 
(AB-EPS08, AB-EPS09, AB-HKN07, and RSZ-LTA). 

STARLIGHT uses the latest HERA data on exclusive J/y production in photon-proton interactions as 
input to calculate the corresponding photon-nucleus cross section. It is worth noting that the HERA 
experimental cross sections have systematic and statistical errors of less than 10% ll2"0"ll2T1l . 

The model by Goncalves and Machado (GM) lfT3l calculates the J/i//-nucleon cross section from the 
Color Dipole model, whereas Cisek, Szczurek, and Schafer (CSS) lfl4l use the essentially equivalent 
k± -factorization approach. The difference of about 25% between the two calculations is due to different 
treatment of the nucleon gluon distribution at low x (gluon saturation), and the way in which it affects 
the dipole-nucleon cross section. 

Calculations by Adeluyi and Bertulani (AB) | fl2l and by Rebyakova, Strikman, and Zhalov (RSZ) 
are based on perturbative QCD. The high mass of the J/y provides a hard scale, and the production can 
to leading order be modelled by two-gluon exchange, with the cross section being proportional to the 
gluon distribution squared. 

The calculations by Rebyakova et al. use a cross section for exclusive J/y photoproduction on a pro- 
ton target, calculated from leading order perturbative QCD within the leading log approximation. The 
calculations use the integrated gluon density distribution in the proton determined by the Durham-PNPI 
group from data on exclusive J/ l/f production at HERA |@). The modification to the nuclear gluon distri- 
bution has been calculated in the Leading Twist Approximation ll29l and is based on using the DGLAP 
evolution equations and the HERA diffractive parton density distributions. 

Adeluyi and Bertulani constrain the nucleon parton distributions to be consistent with data on exclusive 
vector meson production in photon-proton interactions. The photonuclear cross section is then calcu- 
lated using different standard parameterizations of the nuclear gluon distribution functions (EPS08 (AB- 
EPS08), EPS09 (AB-EPS09), HKN07 (AB-HKN07)). For comparison, they also performed calculations 
where the constrained nucleon gluon distribution function is scaled with the number of nucleons without 
shadowing or other nuclear effects (AB-MSTW08). 

In the region of interest here, -3.6< y <-2.6, the sensitivity to shadowing is reduced compared with that 
at mid-rapidity. This is not surprising considering the Bjorken-x ranges probed by this measurement. 
Away from mid-rapidity, there is a two-fold ambiguity in the photon energy and the momentum transfer 
from the nucleus acting as photon target. For example, a. J/y produced at y = 3 corresponds to a photon- 
proton center-of-mass energy of either Wyp =414 GeV or Wyp = 21 GeV. These two energies in turn 
correspond to values of x of about 5 x 10~ 5 and 2 x 10~ 2 , respectively. According to STARLIGHT 
interactions with Wyp = 21 GeV contribute 94% of the cross section, while events with Wyp = 414 GeV 
contribute only 6%. The total da™^/dj at y = 3 is therefore mainly sensitive to the gluon distribution 
around x = 2 x 10~ 2 . 

The measured cross section, dOj°^/dy = 1.00±0.18(stat)^Q 2 g(syst) mb, is compared with the model 
predictions in Fig. 3 a). Fig. 3 b) shows a comparison of the cross section integrated over the range 
-3.6< y <-2.6. The models with largest deviations from the measured value are STARLIGHT and AB- 
MSTW08, which both deviate by about 3 standard deviations if the statistical and systematic errors are 
added in quadrature. Best agreement (within one standard deviation) is seen for the RSZ-LTA, AB- 
EPS09, and AB-EPS08 predictions. A further check can be performed by dividing the rapidity interval 
in two and determining the ratio of the cross sections in each interval. This has the advantage that some 
parts of the systematic errors cancel, and the dominant remaining error is the statistical error. The result is 
R = cr(-3.1 <y < -2.6)/cr(-3.6 <y< -3.1) = 1.36±0.36(stat) ±0.19(syst). The systematic error 
includes the uncertainties in the signal extraction and in the trigger and recontruction efficiency. The 
measured ratio is compared with that from the models in Fig. 3 c). The only models which deviate by 
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more than one standard deviation are AB-MSTW08 and AB-HKN07 (1.7 and 1.5 standard deviations, 
respectively). 

In summary, the first LHC measurement on exclusive photoproduction of J/\j/ in Pb-Pb-collisions at 
■ v /snn = 2.76 TeV has been presented and compared with model calculations. The AB-MSTW08 model, 
which assumes that the forward scattering cross section scales with the number of nucleons squared, 
disagrees with the measurement, both for the value of the cross section and for the ratio of the two 
rapidity intervals, and is strongly disfavoured. STARLIGHT deviates by nearly three standard deviations 
in the cross section and is also disfavoured. Best agreement is found with models which include nuclear 
gluon shadowing consistent with the EPS09 or EPS08 parameterizations (RSZ-LTA, AB-EPS09, and 
AB-EPS08). 



Pb+Pb^ Pb+Pb+Jty ^ = 2.76 TeV 




Fig. 1: Invariant mass distribution for events with exactly two oppositely charged muons satisfying the event 
selection in Table 1, with the mass range extended to 2.2 < M; nv < 4.6 GeV/c 2 . 
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Fig. 2: Dimuon pj distribution for events satisfying the event selection in Table 1, with the px-range extended to 
pj < 0.8 GeV/c. The data points are fitted summing four different Monte Carlo templates: coherent J/y/ production 
(dashed - blue), incoherent J/\j/ production (dotted - red), J/ys from y/ decay (dash-dotted - violet), and yy — >• 
(dash-dotted - green). The solid histogram (black) is the sum. 
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Fig. 3: Measured coherent differential cross section of J/i// photoproduction in ultra-peripheral Pb-Pb collisions 
at ^snn = 2.76 TeV. The error is the quadratic sum of the statistical and systematic errors. The theoretical calcu- 
lations described in the text are also shown. The rapidity distributions are shown in a), b) shows the cross section 
integrated over -3.6 < y < -2.6, and c) shows the ratio of the cross sections in the rapidity intervals -3.1 < y < -2.6 
and -3.6 < y < -3.1. The dashed lines in the lower two plots indicate the three model categories discussed in the 
text. 
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